Non-Cayley-tree model for quasiparticle decay in a quantum dot 
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The decay of a quasiparticle in a confined geometry, resulting from electron-electron interactions, 
has been mapped onto the single-electron problem of diffusion on a Cayley tree by Altshuler et 
al. [Phys. Rev. Lett. 78, 2803 (1997)]. We study an alternative model, that captures the strong 
correlations between the self-energies of different excitations with the same number of quasiparticles. 
The model has a recursion relation for the single-particle density of states that is markedly different 
from the Cayley tree. It remains tractable enough that sufficiently large systems can be studied to 
observe the localization transition in Fock space predicted by Altshuler et al. 
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The lifetime of a quasiparticle in a quantum dot has 
been the subject of recent experimental Q and theoret- 
ical works • Much of the theoretical interest was 
fueled by the striking prediction of Altshuler, Gefen, 
Kamenev, and Levitov || of a critical excitation en- 
ergy below which the lifetime becomes essentially infi- 
nite. This prediction was based on a mapping between 
the decay process of a quasiparticle and the phenomenon 
of Anderson localization on a Cayley tree jl(^, ll) . An in- 
finite lifetime corresponds to the absence of diffusion on 
the lattice in Fock space consisting of n-particle eigen- 
states of the Hamiltonian without interactions. A 
theoretical study of the mapping ||[| predicted a smooth 
transition in the range of excitation energies from Ag x l 2 
to A<? 2 / 3 (with A the single-particle level spacing and g 
the conductance in units of e 2 /h). The thermodynamic 
limit g 1 is essential for the appearance of the transi- 
tion. 

Numerical diagonalizations of a microscopic Hamilto- 
nian H and of the two-body random interaction model 
Q were too far from the thermodynamic limit to observe 
the localization transition. The need for a confirmation 
of the prediction of Altshuler et al. is pressing because 
of a fundamental difference between the decay process in 
Fock space and the diffusion process on a Cayley tree. 
The mapping between the two problems maps different 
'Fn's with the same n onto different sites at the same 
level of the tree. While in the Cayley tree diffusion from 
each of these sites is independent, in the Fock space the 
decay of different 'J'n's is strongly correlated. 

In this paper we consider the model Hamiltonian pro- 
posed by Georgeot and Shepelyansky |0, that permits 
to study these strong correlations in systems that are 
bigger than in Refs. |]||. We find a smooth localization 
transition in Fock space consistent with the predictions 
of Refs. [|]|]J(|. An analytical approximation to our nu- 
merical diagonalizations highlights the origin of the cor- 
relations between the \t„'s. 

The Hamiltonian for spinless fermions is H — Hq + Hi, 
with 



Hi = Vij,kl4 c k c i c j- W 

i<j,k<l 



The non-interacting part Hq contains the single-particle 
levels £j in a disordered quantum dot. We count the lev- 
els from the Fermi level, meaning that the ground state 
of Hq has occupied levels for j < and empty levels 
for j > 0. We assume that an energy level £j is uni- 
formly distributed in the interval [(j — |)A, (j + |)A]. 
This yields a linear level repulsion, consistent with time- 
reversal symmetry. The basis of Hq consists of states 
that have m electron excitations (occupied levels with 
j > 0) and n hole excitations (empty levels with j < 0). 
The two-body interaction Hi couples them to states that 
differ by at most two electron-hole pairs. 

We assume that Vy,y = 0. (These diagonal terms can 
be incorporated into Hq in a mean- field approximation.) 
For the off-diagonal matrix elements we adopt the layer 
model of Ref . [Q , which is based on the following obser- 
vation. The interaction strength V is related to A and 
g by HP V = A/g. Since V « A for g > 1, only 
eigenstates of Hq within an energy layer of width A are 
strongly coupled by the interaction. The layer model ex- 
ploits this in a clever way by setting VijM = unless 
i, j, fc, I are four distinct indices with i + j = k + I. The 
non-zero 14j,fcz are chosen to be independent real random 
variables, subject to the restriction Vij,hi = Vki,ij im- 
posed by the hermiticity of the Hamiltonian. We also set 
Vji,ki — —Vij,kl — Vij,ik- The distribution of each matrix 
element is taken to be a Gaussian with zero mean and 
variance V 2 . 

One advantage of the layer model is that the ground 
state |FS) of Ho (the Fermi sea) remains an eigenstate 
of Hq + Hi. We assume that it remains the ground 
state. A second advantage is that the effective dimen- 
sion of the Hilbert space is greatly reduced. The num- 
ber of states into which an electron excitation cj|FS) 
of energy Ej decays is equal to the number V(j) ~ 
(4j\/3) -1 exp(7r^/2j/3) of partitions of j, for a suffi- 
ciently large number of electrons in the quantum dot [jl4| . 
This grows much more slowly with j than in the conven- 
tional two-body random interaction model [jl5 16 1, used 
in previous work P,|9|,|l7t . While the layer model renders 
the problem tractable, it preserves the strong correlations 
mentioned in the introduction, as we will discuss shortly. 
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The decay of the quasiparticle state cj|FS) is described 
by the Green function 

G j (E) = (FS\c j (E 
= [E~ e, - 



+ E FS -H) V 



]|FS) 



(2) 



where E-ps is the energy of the Fermi sea: i?|FS) = 
F |FS) = £fs|FS). The second equality in Eq. (§) 
defines the self-energy Ej(E). The quantity of physi- 
cal interest (measured by means of a tunneling probe in 
Ref. is the single-particle density of states Pj{E) — 
J2 a S(E + Efs — -E Q )|(a|cj|FS)| 2 , where the sum over a 
runs over all eigenstates \a) of H, with eigenvalues E a . 
It is related to the imaginary part of the Green function 
by 



(E) = V^ kl \E ~ it - S k - e t - E,(£? - S k - e,) 



A7 



E fc (£7-ej -Ei) -^i{E-S l -S k ) 



(7) 



where the sum runs over the indices k, I with < k < I 
and i = k + I — j < 0. 

This recursion relation in Fock space can be compared 
with the recursion relation for the Cayley tree [jl0| , which 
has the form 



e k -E k (E)Y 



(8) 



pj(E) = lim Im GAE + irj). 
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The ensemble average Pj(E) is not sensitive to the der- 
ealization transition. For that reason, we will also study 
the inverse participation ratio Pj(E) — J2 a ^C® + -^fs — 
E a )\ (a\c!j |FS) | 4 , related to the Green function by 

P j (E) = -limr ] \G j (E + ir ) )\ 2 . (4) 
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The dimensionless ensemble-averaged quantity 

Pj = Pi(zj)/pJ&) (5) 

increases from to 1 on going from extended to localized 
states. 

We have computed the Green function numerically us- 
ing an iterative Lanczos method. The largest system we 
could study in this way has j — 25, corresponding to 
a basis of r P(25) — 1958 states. Before presenting the 
results of this exact diagonalization, we discuss a cer- 
tain decoupling approximation that has the advantage of 
showing explicitly how the decay of the quasiparticle is 
different from the diffusion on a Cayley tree. 

The problem of the diffusion on a Cayley tree can be 
solved exactly because the self-energy satisfies a closed 
recursion relation pC|JTl[ . Such a recursion relation ex- 
ists because the Cayley tree has no loops. The lattice in 
Fock space generated by the quasiparticle decay process 
H does have loops, but we believe that these do not play 
an essential role and we will ignore them. The decoupling 
approximation consists in writing the self-energy E 4 fej (E) 
of a three-particle excitation as the sum of single-particle 
self-energies: 



Z m (E) = Zi(E- e k - et) + E fc (E - 
+ Y, l {E-e i -e k ). 



£1 - £i) 



(0) 



Here is the excitation energy, defined as Si = £j for an 
electron (i > 0) and Si — —Ei for a hole (i < 0). With 
this approximation, the self-energy satisfies the recursion 
relation 



Here the sum runs over all sites k (energy e k ) of the next 
level of the tree that are connected to j, with hopping 
matrix elements Vj ik . We notice two differences between 
Eqs. (Q) and (g). The first is that the recursion relation 
on the Cayley tree conserves energy, while the recursion 
relation in Fock space does not. Another way of say- 
ing this is that Eq. (|J) is a recursion relation between 
numbers at one fixed E, while Eq. (|7j) is a relation 
between functions T,j(E). The second difference is that 
the number of self-energies coupled by repeated applica- 
tions of the recursion relation in the Cayley tree grows 
exponentially (limited only by the size of the lattice), 
while in Fock space this number remains fixed at the 
number 2j of single-particle levels coupled to the excita- 
tion ct|FS) by the interaction. Since 2j is exponentially 
smaller than the size V{j) of the lattice in Fock space, 
this is an enormous difference with the Cayley tree. We 
are able to make such a precise statement because of the 
simplifications inherent to the layer model. However, we 
believe that the strong correlations between excitations 
with the same number of quasiparticlcs implied by Eq. 
(0) are present as well in the full problem of quasiparticle 
decay — although we can not write down such a simple 
recursion relation for the full problem. 

We have calculated the average single-particle density 
of states pj (E) and the inverse participation ratio Pj by 
exact diagonalization for j up to 25, as a function of 
the dimensionless conductance g = A/V. (We have also 
computed the same quantities by numerically solving the 
recursion relation (|7]), and find good agreement.) The 
results for the average density of states collapse approxi- 
mately onto the same curve (see Fig. [l]) , once the energies 
are rescaled by T = ^TrA(j / g) 2 . This expression for T is 
the decay rate following from the golden rule of pertur- 
bation theory Q , assuming an energy-independent three- 
particle density of states (equal to |j 2 /A in the layer 
model). The small deviations from a Lorentzian (dotted 
curve in Fig. |l|) are an artefact of the layer model. (They 
disappear if the restriction i + j = k + / on the matrix 
elements Vi^u of the interaction is removed.) 
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Fig. [j] — Average single-particle density of states pj (E) , 
rescaled by T = ^nA(j / g) 2 , for j = 25. The solid and 
dashed curves are computed by exact diagonalization of 
the layer model for g = 55 and g = 300, respectively. 
Averages are taken over 7500 realisations of the random 
Hamiltonian. The dotted curve is a Lorentzian of unit 
area and width. 

As expected, there is no indication in the average den- 
sity of states pj of a localization transition. The density 
of states pj for a single realisation of H is shown in Fig. 
||. The difference between pj for small and large values 
of the ratio T/V is striking, and in qualitative agreement 
with the prediction of Altshuler et al. ||] : Sharp isolated 
peaks in the localized regime (top panel), in contrast to 
a single broad peak in the delocalized regime (bottom 
panel) . 

To study the localization transition we calculate the 
inverse participation ratio Pj , defined in Eq. (|s|) . Follow- 
ing Ref. |)|, we compare with the prediction of a totally 
delocalized situation ("golden rule"). The golden rule 
prediction is Pj ~ min(l,(5/r), where 5 is the mean en- 
ergy separation of the eigenstates |a) of H . In the layer 
model, 8 ~ A/V(j). Since S/T oc g 2 , the golden rule 
predicts a quadratic increase of Pj with increasing g , 
until Pj saturates at a value of order unity. A faster 
than quadratic increase is a signature of localization. We 
show in Fig. [| a double-logarithmic plot of Pj versus g 
for j = 15, 20, and 25. The straight lines of slope 2 show 
the quadratic increase predicted by the golden rule. The 
largest system considered (j = 25, squares) has unam- 
biguously a region of faster than quadratic increase of 
Pj, starting at g w 100 and persisting until saturation is 
erached at g ss 500. In contrast, the smallest system con- 
sidered (j — 15, triangles) follows the golden rule predic- 
tion until it saturates at g sa 60. This system is clearly 
too small to show the transition to a localized regime. 
The largest system studied in Ref. |jj had j ~ 15, and 
indeed no deviations from the golden rule were found in 
that paper. 
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Fig. H — Single-particle density of states pj of an individ- 
ual member of the ensemble of quantum dots, computed 
by exact diagonalization for j = 25 and two values of 
g : 300 (upper panel) and 55 (lower panel). The two 
results are qualitatively different, although the ensemble 
averages are essentially the same (see Fig. [l]). 
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Fig. [I — Inverse participation ratio as a function of di- 
mensionless conductance, for j = 15 (triangles), j = 20 
(circles), and j = 25 (squares). The straight lines of 
slope 2 on the log-log scale show the quadratic increase 
of Pj with g predicted by the golden rule. A faster than 
quadratic increase indicates a transition to the localized 
regime. Statistical error bars have the size of the mark- 
ers. 
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In conclusion, we have studied a model for quasipar- 
ticle decay in a quantum dot that preserves the strong 
correlations omitted in the Cayley tree model, yet re- 
mains tractable enough that large excitation energies are 
accessible. Our largest system demonstrates a transition 
to a localized regime that had remained elusive in previ- 
ous studies on smaller systems. 
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